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Abstract 

The properties of the low-lying, collective 2f states in neutron-rich oxygen 
isotopes are investigated in the framework of self-consistent microscopic mod- 
els with effective Skyrme interactions. In RPA the excitation energies E 2 + 
can be well described but the transition probabilities are much too small as 
compared to experiment. Pairing correlations are then accounted for by per- 
forming quasiparticle RPA calculations. This improves considerably the pre- 
dictions of B(E2) values and it enables one to calculate more reliably the ratios 
M n /M p of neutron-to-proton transition amplitudes. A satisfactory agreement 
with the existing experimental values of M n /M p is obtained. 

PACS numbers: 23.20. Js, 21.60.Jz, 21.10.Re 
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The prospects of nuclear physics studies with nuclei far from stability open up a wide 
range of possibilities for refining our understanding of nuclear properties in terms of mi- 
croscopic descriptions and effective nucleon-nucleon interactions. One of the important 
aspects is the ability to disentangle neutron and proton contributions to collective transi- 
tions between low-lying states and the ground state. Experimentally, this can be done in a 
phenomenological way [|J by combining the information obtained in measurements involving 
various hadronic probes and purely electromagnetic probes. For instance, numerous exper- 
imental studies have been performed on 18 using different hadronic probes like proton 0], 
neutron ||, or pion scattering. More recently, proton scattering on 20 O |J yielded infor- 
mation on the first 2 + state in this unstable neutron-rich oxygen isotope. On the other hand, 
studies involving only electromagnetic properties such as electron scattering @, Coulomb 
excitation or lifetime measurement [0] have also been done for these nuclei. While the exci- 
tation processes of purely electromagnetic nature are sensitive only to the protons and give 
access to the proton transition amplitudes and transition densities, the hadronic processes 
are sensitive to both proton and neutron transition densities. Therefore, it is possible by 
a combined analysis of the data from electromagnetic and hadronic processes to determine 
experimentally for a given excited state the transition amplitudes M v and M n corresponding 
to protons and neutrons, respectively HJ. Proton scattering experiments yielding M n /M p 
values have been recently performed on neutron-rich sulfur and oxygen isotopes |6[,|9|, |l0|] . 

In this work, we present microscopic calculations of low-lying 2 + states in neutron-rich 
oxygen isotopes. These calculations are based on effective Skyrme interactions and they are 
performed in the framework of the random phase approximation (RPA) and the quasiparticle 
random phase approximation (QRPA). In microscopic models the properties of the states de- 
pend on two main inputs, the single-particle spectrum and the residual two-body interaction. 
In the present approach these two features are linked since the same effective interaction 
determines the Hartree-Fock (HF) single-particle spectrum and the residual particle-hole 
interaction. This approach has proved to be an efficient mean to predict properties of collec- 
tive excitations like giant resonances |nj and it has also been used for calculating low-lying 
collective states in closed-shell nuclei [ |12|1 . 

In unstable nuclei we usually don't deal with closed-shell or closed-subshell systems and 
therefore, the HF and RPA calculations must be done with additional approximations. The 
HF calculations are carried out assuming spherical symmetry and using the standard filling 
approximation with equal occupation numbers for all (jm)-substates of the partially filled j- 
subshell. The RPA calculations are then be carried out taking into account these occupation 
numbers. However, pairing correlations can be important in such nuclei and they must be 
taken into account. Their effects will be described by HFBCS calculations for the ground 
states and by QRPA calculations for the excitation spectra. 

The HF and HFBCS method in spherical nuclei with Skyrme interactions is well-known 
i3| , |T4|| . For the pairing interaction we simply choose a constant gap given by ||15|| : 



A = 12.A~sM.eV . (1) 

In a more realistic treatment of the pairing, the gap would depend on the single-particle 
state considered and it would tend to zero when the subshell is far from the Fermi level. 
Thus, in the constant gap approximation it is necessary to introduce a cut-off in the single- 
particle space. Above this cutoff subshells don't participate to the pairing effect. In the case 
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of oxygen isotopes, we choose the BCS subspace to include the Is, lp and 2s — Id major 
shells. 

The results of the HF and HFBCS calculations performed for the nuclei 18 > 20 > 22 O using 
typical Skyrme interactions are summarized in Table |. For all the interactions, the binding 
energies per particle decrease with increasing neutron number but SGII has a different behav- 
ior as compared to the other interactions and it predicts larger values of B/A. The pairing 
correlations decrease B/A by about 4% in all cases. The neutron radii increase substantially 
from 18 O to 22 O and the proton radii also increase slightly as a result of neutron-proton 
attraction. The effect of neutron pairing correlations is to redistribute neutron densities to 
the tail region and therefore, this leads to a small increase in r n (and also in r p for the reason 
mentioned above). 

The HF-RPA model with Skyrme effective forces is also well-known. We only mention 
that in this work we solve the RPA equations in configuration space, choosing the particle- 
hole space so as to exhaust the energy-weighted sum rules. The continuous part of the 
single-particle spectrum is discretized by diagonalizing the HF hamiltonian on a harmonic 



oscillator basis ||17|| . To generalize the HF-RPA to the QRPA model we follow the standard 
procedure Denoting by a^, a a the creation and annihilation operators of a particle in 
a HF state a = (j a , m a ) and by c^, c a the corresponding operators for a quasiparticle state, 
we have: 



C ] a m a U a a j a m a V a 



(-l) ja+ma a ja - ma , (2) 



where the BCS amplitudes u a , v a satisfy the normalization condition: 

ul + v 2 a = l. (3) 

These amplitudes are determined by solving the BCS equations. One can then build the 
two-quasiparticle creation operators in an angular momentum coupled scheme: 

CyjM) = (1 + M~ 1/2 E Uajpm^l JM)4 c t . (4) 

m a mp 

In QRPA the nuclear excitations correspond to phonon operators which are linear combina- 
tions of two-quasiparticle creation and annihilation operators: 

Q V \JM) = £ X^(J)C^(JM) + {-l) M Y^{J)C aP {J - M) . (5) 

a>/3 

Making use of the condition that the QRPA ground state |0) is a vacuum of phonon: 

Q v (JM)\0) = 0, (6) 

one can then derive the QRPA equations whose solutions yield the excitation energies E v 
and amplitudes X^s, Y£p of the excited states. 

An important quantity that characterizes a given state v = (E u , LJ) is its transition 
density: 

<^(r)^(H£<H r - r i)|o> , (7) 
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and a similar definition of the neutron (proton) transition density bp v n (Spp) with the sum- 
mation in eq.(0) restricted to neutrons (protons). In QRPA the radial part of the transition 
density is: 

8p v {r) = W(r) < flinolla > {X^(J) - Y^(J)}{u a v p + {-l) J v a u p } , (8) 

a>/3 

where <£> a (r) is the radial part of the wavefunction of the quasiparticle state a. As an example 
the QRPA neutron and proton transition densities of the first 2 + state in 20 O calculated with 
the interaction SGII are shown in Fig.l. The neutron transition density is shifted outwards 
as compared to the proton transition density due to the presence of a neutron skin. Clearly, 
the two transition densities do not scale like N/Z as it is sometimes assumed and this has 
quantitative consequences as we shall see below. 

The neutron and proton matrix elements M =< z/|r L Y^ |0 > of a multipole operator are 
obtained by integrating the corresponding transition densities over r : 

M n , P = J 5p» iP (r)r L+2 dr , (9) 
and the reduced electric multipole transition probabilities are calculated as 

B{EL) n , p = |M„, P | 2 . (10) 

We have calculated the J n = 2+ states in 18 > 20 > 22 O using RPA and QRPA with SIII, SGII 
and SLy4 interactions. The energies and B(E2) P values of the first 2 + states are shown 
in Fig. 2 together with the existing experimental values in 18 O and 20 O. The data come 
from experiments involving electromagnetic processes such as Coulomb excitation or lifetime 
measurement ||. For the E 2 + energies, standard RPA reproduces very well experimental 
values especially for SGII and SLy4. On the other hand, QRPA deteriorates this agreement 
and it predicts the first 2 + states at somewhat higher energies. The theoretical prediction of 
the energies of low-lying states in models based on a HF or HFBCS mean field is a delicate 
task because these energies are sensitive to the spin-orbit part of the mean field while the 
spin-orbit component of the two-body effective interaction is not so well determined. For 
the three interactions used here there is a clear difference in the QRPA E 2 + energies between 
SGII and the other interactions. In the case of B(E2) P values RPA predicts too small values 
with the three interactions. However, the agreement becomes very good for QRPA both in 
18 and 20 O, especially for SGII and SLy4 interactions. This indicates that pairing effects 
are important for transition probabilities of the first 2 + state in these nuclei. 

In Fig. 3 are shown the ratios M n /M p calculated with the three interactions within RPA 
and QRPA. On the same figure are displayed the experimental values taken from ref ||. 
It can be seen that the RPA results are somewhat larger than those of QRPA due to 
the very small B(E 2 ) P values obtained in RPA. In comparison with the data, the QRPA 
results are in good agreement in 20 O whereas they are slightly too large in 18 0. The three 
interactions predict different A dependence in these oxygen isotopes for the QRPA ratios. 
It would be interesting to obtain the experimental value of M n /M p in 22 as well as the 
corresponding B(E2) transition probability. If one assumes the quadrupole excitation to 
be purely isoscalar the ratio M n /M p would be equal to N/Z. Taking as a guideline the 
QRPA results calculated with SGII one sees that the ratio of neutron-to-proton transition 
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amplitudes is about (1.8 - 2.0)N/Z, thus indicating that the low-lying 2 + state has an 
important isovector component. This is not surprising since in these neutron-rich nuclei 
there are neutron particle-hole configurations at low energy which have no counterpart on the 
proton side and therefore, these neutron configurations necessarily introduce both isoscalar 
and isovector type of excitations. 

In summary, we have investigated the properties of the low-lying, collective 2f states in 
neutron-rich oxygen isotopes in the framework of self-consistent microscopic models. Within 
the RPA model the excitation energy E 2 + can be well described but the transition probabili- 
ties are much too small as compared to experiment. In these open subshell nuclei the pairing 
correlations can be important and therefore, we have extended the previous model to the 
HFBCS approximation at the mean field level and we have performed QRPA calculations 
for the excited states. The quasiparticle microscopic description improves considerably the 
predictions of B(E2) values and it enables us to calculate more reliably the ratios M n /M p of 
neutron-to-proton transition amplitudes. These calculated values differ noticeably from the 
naive N/Z estimate and they are in satisfactory agreement with experiment. However, the 
QRPA overestimates the E 2 + energies. The sensitivity of positive parity low-lying states to 
the effective interaction should give one a handle on some specific components of the force, 
for instance the two body spin-orbit part. 

Further proton scattering results on 18 and 20 O will be available soon [HJ, yielding 
energies and M n /M p ratio for the first 2 + and 3~ states. It would also be useful to perform 
experiments on more neutron- rich oxygen isotopes to establish firmly the trend of M n /M p 
function of N/ Z. 

We would like to thank G. Colo, T. Suomijarvi and C. Volpe for useful discussions. 
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Figure captions 



Figure 1. Neutron and proton transition densities of the first 2 + state in 20 O, calculated 
in QRPA with interaction SGII. 

Figure 2. Energies and B(E2) p values of the first 2 + states in oxygen isotopes. Open and 
black symbols correspond to RPA and QRPA calculations, respectively. Three effective in- 
teractions are used: SIII |TJ[ (circles), SGII |Tl| (triangles), SLy4 [TfJ (stars). Experimental 
values are shown as crosses, with error bars for B(E2) P . 

Figure 3. The M n /M p ratios in oxygen isotopes. The notations are the same as in Fig. 

2. 
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TABLES 







18Q 


20Q 


22Q 


r n (fm) 


r p (fm) 


B/A (MeV) 


r n (fm) 


ip (fm) 


B/A (MeV) 


T n (fm) 


t p (fm) 


B/A (MeV) 


SIII 


HF 


2.77 


2.65 


7.74 


2.89 


2.67 


7.58 


2.97 


2.69 


7.48 


BCS 


2.79 


2.67 


7.43 


2.91 


2.69 


7.28 


3.03 


2.71 


7.16 


SGII 


HF 


2.76 


2.64 


8.30 


2.87 


2.65 


8.15 


2.95 


2.66 


8.07 


BCS 


2.80 


2.66 


7.95 


2.92 


2.68 


7.79 


3.02 


2.69 


7.69 


SLy4 


HF 


2.83 


2.69 


7.74 


2.94 


2.70 


7.57 


3.02 


2.71 


7.46 


BCS 


2.85 


2.72 


7.42 


2.98 


2.73 


7.25 


3.08 


2.74 


7.13 



TABLE I. 

Neutron and proton r.m.s. radii, and binding energy per particle in the nuclei 18 0, 20 O and 
22 calculated with interactions SIII Q, SGII Q and SLy4 @. The rows labeled HF 
and BCS show Hartree-Fock and Hartree-Fock-BCS results, respectively. 
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